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INVESTIGATIONS  IN  SUPPORT  OF 
HIGH  ENERGY  LASER  TECHNOLOGY 

lo  0  INTRODUCTION 

During  the  past  several  years,  a  series  of  analytical  investiga¬ 
tions  into  the  physics  and  chemistry  of  various  laser  media  has  been 
performed.  Early  in  the  program,  line-width  measurements  of  HF 
broadening  by  various  gases  were  analyzed,  and  the  results  were  used 
to  do  an  extensive  set  of  calculations  of  both  HF  and  DF  line  widths 
at  chemical  laser  cavity  conditions  for  gain-profile  modeling.  This 
study  of  the  broadening  of  laser  molecular  line  profiles  led  to  an  in¬ 
vestigation  of  inhomogeneous  broadening  and  cross  relaxation  of  the 
laser  gain  profile.  Recent  results  of  this  study  are  in  this  report. 

Also  included  here  are  calculations  of  the  vibration  deactivation  of 
HCl  by  H  atoms.  This  study  is  significant  for  modeling  the  promis¬ 
ing  non-reacting  mixing  laser  which  would  operate  on  the  HCl  vibra¬ 
tion  transitions  in  the  infrared. 

Investigations  of  potential  visible  laser  molecules  have  been  a 
continuing  part  of  the  present  program.  The  visible  laser  studies 
have  been  fundamental  in  nature  and  have  mainly  addressed  the  alkali- 
noble  gas  laser  media.  Results  of  these  calculations  are  included  in 
the  present  volume  along  with  some  applications  of  the  inhomogeneous 
broadening  formulation  to  visible  laser  systems. 
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2.  0  H  ATOM  DEACTIVATION  OF  HCl  VIBRATION 

Lasing  from  the  vibrationally  excited  states  of  HCl  has  been 
achieved  by  a  variety  of  excitation  mechanisms.  These  mechanisms 
produce  vibrationally  hot  and  rotationally  cold  HCl;  thereby  achieving 
lasing  in  the  P  branch  of  the  vibrational  transitions. 

The  maintenance  of  this  disparity  between  the  vibrational  and 
rotational  temperatures  is  paramount  to  device  performance.  Cooling 
processes  inherent  to  the  laser  medium  or  to  the  excitation  mechanism 
limit  achievable  device  performance.  One  formidable  coolant  present 
in  HCl  laser  media  is  the  H  atom.  This  specie  is  always  produced 
by  the  rapid  dissociation  of  HCl  arising  from  its  elevated  vibrational 
temperature.  Several  excitation  mechanisms  produce  additional  H 
atoms,  e.  g.,  the  chain  reactions  operative  in  the  chemical  laser 
devices  require  the  presence  of  H  atoms.  In  the  nonreacting  mixing 
laser  [1],  H  atoms  are  produced  in  the  vibrational  heating  of  the  H2 
pumping  gas.  Hydrogen  atoms  are  generally  present  in  HCl  laser  media. 
Models  of  proposed  high-power  HCl  lasers  show  the  effect  of  H  atoms 
to  be  a  dramatic  reduction  in  performance;  therefore,  their  cooling 
effects  must  be  considered  in  the  engineering  of  high  performance  HCl 
laser  devices. 

Available  data  is  insufficient  for  the  modelling  of  high  energy 
devices.  The  room  temperature  rate  of  the  process 

H  +  HCl  (v=l)  — -^H  +  HCl  (v=0)  (1) 

has  been  measured  [2].  Theoretical  attempts  [3]  to  extrapolate  this  datum 
to  other  conditions  have  been  unreliable,  and  only  qualitative  empirical 
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extrapolations  have  been  employed  in  device  modelling.  The  present 
study  provides  a  reliable  theoretical  determination  of  the  temperature 
dependence  of  the  measured  process.  In  addition,  the  detailed  rates 
of  deactivation  of  the  higher  states  v=2  and  3  are  related  to  the  v=l 
to  v=0  process. 

2. 1.  Potential  Energy  Surfaces 

As  pointed  out  above,  there  is  a  dearth  of  data  on  the  colli¬ 
sions  of  H  atoms  with  HCl.  (Extrapolation  from  a  single  data  point  is 
hardly  justified. )  Fortunately,  considerable  data  is  available  on  the 
related  reactive  process 

Hg  +  Cl  — ^  HCl  +  H.  (2) 

These  data  have  been  analyzed  to  extract  the  potential  operative  during 
reactive  encounters  [4],  a  related  aspect  of  the  potential  controlling 
collisions  of  H  atoms  with  HCl  molecules.  From  these  data,  the 
potential  controlling  collisions  of  the  atom  with  the  H  end  of  the  mole¬ 
cule  is  known.  Attack  on  the  opposite  end  of  the  molecule  samples 
regions  of  the  potential  which  play  no  role  in  the  reaction;  hence, 
those  segments  of  the  potential  are  unknown. 

In  the  absence  of  adequate  empirical  data  on  this  region  of 

the  surface,  calculation  of  the  potential  is  a  viable  alternative  for 

the  determination  of  the  potential.  The  Diatomics-in-Molecules  [5] 

potential  for  this  region  was  found  to  exhibit  an  11  kcal/mole  barrier 

to  the  H  atom  exchange  process  in  the  linear  symmetric  geometry  at 

o 

an  H-Cl  separation  of  1.  46A.  Similar  calculations  [6]  show  a  16  kcal/ 
mole  barrier  in  the  F  analogue  of  this  system.  In  contrast,  accurate 
ab  initio  calculations  [7]  indicate  the  barrier  there  is  iii  excess  of  40 
kcal/mole.  Thus,  the  DIM  barrier  height  is  suspect,  and  the  actual 
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barrier  may  be  considerably  higher  than  the  11  kcal/mole  DIM  pre¬ 
diction.  Scaling  the  ab  initio  barrier  in  H  F  H  using  the  DIM  formalism 
suggests  a  28  kcal/mole  barrier  is  present  in  the  H- Cl -H  system. 

The  LEPS  potential  function  used  in  the  scattering  calculation  exhibits 
an  8  kcal/mole  well  in  this  region  [8]  and  is  clearly  erroneous  there. 
Collisions  striking  the  Cl  end  of  the  molecule  will  be  predicted  to 
undergo  direct  H  atom  exchange  rather  than  simple  elastic  collision. 
(Mass  factors  make  inelastic  processes  improbable. )  For  the  analysis 
these  collisions  have  simply  been  neglected  to  simulate  the  effect  of 
the  missing  barrier;  the  inelastic  processes  of  interest  occur  in  the 
known  region  of  the  potential.  Therefore,  the  empirical  potential 
is  adequate. 

2.  2.  Kinetics  Calculations 

Classical  scattering  calculations  on  this  potential  surface 
have  been  used  to  estimate  the  kinetics  of  H  atom  deactivation  of  HCl 
vibration.  These  calculations  have  been  performed  using  an  exten¬ 
sively  modified  version  of  Polanyi’s  REACT2  code  [9].  These  revisions 
included  recoding  of  the  numerical  integration  sections  to  improve 
the  efficiency  of  trajectory  tracking,  changes  in  the  Monte-Carlo 
selection  of  representative  collisions,  and  insertion  of  an  analysis 
routine  specific  to  the  present  problem. 

Analysis  of  the  calculations  yields  rate  data  pertinent  to 
modelling  of  HCl  laser  cavity  chemistry.  The  thermal  rate  of  reac¬ 
tive  consumption  of  HCl  by  H 

H  +  HCl - Hg  +  Cl  (3) 

is  found  to  be  slow  (1. 0  x  10"^^  cm^/sec  at  lOOO^K);  however,  this 
rate  rises  exponentially  with  HCl  vibrational  state  (see  Figure  1). 
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The  computed  room  temperature  rate  constant  of  the  H  atom 
deactivation  of  the  v=l  state  of  HCl,  (2. 4  ±  .  4)  x  cm  /sec  is 

in  fair  agreement  with  the  measured  rate  constant,  (6.  5  —  2.  2)  x  10 
cm^/sec.  Elimination  of  the  atom  exchange  deactivation  process  re¬ 
duces  the  rate  constant  from  Wilkin’s  [3]  value  of  5.0  x  10  cm  /secj 
the  remaining  process  is  clearly  the  observed  inelastic  process. 

The  computed  rate  constant  is  overestimated  because  of 
inaccuracies  in  the  empirical  potential  employed.  The  potential 
governing  perpendicular  H  atom  attack  is  not  well  specified  by  the 
data  employed  in  determining  this  potential;  therefore,  sizeable  errors 
may  occur  in  that  region.  By  continuity,  the  underestimate  of  the 
unknown  potential  for  attack  on  the  Cl  end  of  the  molecule  cited  above 
implies  an  underestimate  for  the  poorly  known  region  here.  For  this 
reason,  the  solid  angle  in  which  H  atom  attack  is  effective  should  be 
overestimated.  Concommitantly,  the  rate  of  the  deactivation  reaction 
should  be  overestimated,  but  the  fundamental  physics  of  the  process 
described  is  correct.  Therefore,  the  trends  predicted  in  this  study- 
temperature  and  state  dependences  —  are  expected  to  be  much  more 
accurate  than  the  magnitudes  of  the  computed  effects. 

2. 3.  Results 

The  techniques  described  above  have  been  applied  to  the 
rate  of  H  atom  collisional  cooling  of  vibrationally  excited  HCl 

H  +  HCl  (v) - »-H  +  HCl  (v’ )  (v'  <  v).  (4) 

The  temperature  dependence  of  the  1  0  process  has  been  mapped 

from  300®K  to  1000°K,  and  the  vibrational  state  dependence  has  been 
computed  at  500^K. 
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The  rates  of  the  1 — ►  0  process  given  in  Table  I  suggest  a 
negative  activation  energy,  -316  cal/mole.  Within  the  accuracy  of 
the  present  calculations,  this  rate  constant 

Kj  q(T)  =  AT^  (300°<T<1000°)  (5) 

is  indistinguishable  from  the  activationless  process.  Therefore,  the 
assumption  of  that  form  introduces  no  significant  error  in  modelling 
laser  cavity  chemistry. 

The  relative  rates  of  the  deactivation  of  the  higher  states  are 
compared  in  Figure  1  and  Table  II.  These  results  clearly  show  the 
inapplicability  of  the  cannonical  assumption  that  multiquantum  colli- 
sional  cooling  is  negligible.  To  the  contrary,  for  this  system  the 
higher  Av  processes  are  faster  than  those  for  lower  Av's.  The  rates 
of  the  detailed  deactivation  processes  appear  to  decrease  linearly 
with  initial  vibrational  state  up  to  v=3.  (Beyond  this  value,  contribu¬ 
tions  from  the  atom  exchange  process  should  alter  the  shapes  of  the 
curves  given. ) 

The  nature  of  the  deactivation  process  is  described  by  the 
mean  change  in  translational  energy  associated  with  the  deactivation. 
Figure  2  shows  the  v— ^0  processes  is  essentially  pure  v — ^T  transfer. 
Transitions  into  the  higher  vibrational  states  leave  increasing  amounts 
of  energy  in  diatomic  rotation,  but  the  process  remains  predominantly 
v-^T  except  for  the  3-^2  process.  Even  in  that  case,  half  the  energy 
is  deposited  in  the  translational  mode. 

The  loss  of  HCl  vibrational  excitation  to  H  atom  collisions 
involves  not  only  the  deactivation  processes  described  above,  but  also 
the  reaction  to  form  H2 

H  +  HCl — — +  Cl.  (6) 
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Table  1.  Rate  constants  for  H  atom 
deactivation  of  HCl  (v=l) 
between  300^  and  lOOO^K. 


T  (°K) 

Q(10"^^cm^/sec) 

300 

2.35  ±  .39 

500 

2.38  ±  .16 

750 

2.  00  ±  .  46 

1000 

2.97  ±  .12 

Table  II.  Relative  rate  constants  for 
H  atom  deactivation  of 
HCl  at  BOO^K 


V 

V’ 

1 

0 

1.00 

2 

1 

0.  62 

2 

0 

0.  80 

3 

2 

0.  23 

3 

1 

0.48 

3 

0 

0.  67 
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Although  slightly  exothermic,  this  reaction  is  negligible  under  thermal 
conditions  because  it  is  inhibited  by  mass  factors.  The  results  in 
Figure  1  show  the  reaction  rate  rises  exponentially  with  the  HCl 
vibrational  state  to  a  value  comparable  to  the  deactivation  processes 
for  v=3.  This  leads  to  a  systematic  depopulation  of  the  higher  states, 
and  therefore  has  a  deleterious  effect  on  laser  performance. 

State  population  is  reduced  by  H  atom  deactivation  to  all  lower 
states  and  by  chemical  reaction  with  those  H  atoms.  The  total  rate  of 
depopulation  of  the  individual  vibrational  states  is  shown  in  Figure  1. 
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3.  0  THE  EFFECT  OF  INHOMOGENEOUS  BROADENING  AND 
CROSS  RELAXATION  ON  HEL  SYSTEM  PERFORMANCE 

A  theoretical  framework  was  provided  in  a  previous  report  [10], 
to  assess  the  importance  of  inhomogeneous  broadening  and  cross  re¬ 
laxation  on  high  energy  laser  system  performance.  In  the  present 
report  the  results  of  the  previous  work  are  extended,  and  the  theory 
is  applied  to  some  additional  laser  systems  of  current  interest. 

3. 1  Theory  Review 


The  basic  parameter  f^^^  is  used  to  assess  adequacy  of  the 
homogeneous  line  shape  assumption  for  a  laser  oscillator  [11, 12]. 
This  parameter  is  defined  as  follows: 


f 

osc 


(7) 


When  f^g^  is  found  to  be  greater  than  one,  it  may  be  safely  assumed 
that  the  cross -relaxation  rate  and  homogeneous  line  width  combination 
is  large  enough  so  that  the  line  may  be  assumed  to  be  homogeneously 
broadened.  If  however  f^^^  is  less  than  one,  the  homogeneous  line 
assumption  may  not  be  valid.  In  the  later  case  a  more  careful  con¬ 
sideration  of  the  particular  laser  system  would  be  warranted  to  de¬ 
termine  if  a  hole-burning  type  of  gain  saturation  is  limiting  the  output 
power  of  the  laser. 

Formulas  for  the  various  parameters  occurring  in  equation  (7) 
are  repeated  below  because  there  are  some  typographical  errors  in 
the  original  report  in  the  definition  of  these  quantities.  The  cross - 
relaxation  rate  is  assumed  equal  to  the  molecular  collision  rate  and 
is  given  by 

F  =  4.6924*  (sec'^)  (8) 

n 


Here  a  is  the  kinetic  theory  diameter  in  A,  p  is  the  total  cavity  pres¬ 
sure  in  atmospheres,  M  is  the  reduced  mass  of  the  active  molecule 
and  perturber,  and  T  is  the  Kelvin  temperature.  S  is  the  stimulated 
emission  rate  per  upper  state  molecule.  The  Lorentzian  half -width 
is  given  by 


0.1766  cr’‘ 


(cm"^)  (9) 


where  in  the  above  a'  is  the  effective  collision  diameter  for  line 
broadening.  The  standard  form  of  the  Doppler  half  width  is  assumed 
(this  equation  was  in  error  in  the  previous  report); 


D 


3.58  X  V 


0 


(cm”^)  (10) 


In  this  expression,  is  the  frequency  of  the  transition  in  wave- 
-1  ^ 

number  (cm  ),  and  m  is  the  mass  of  the  active  molecule  in  atomic 
mass  units.  Substitution  of  the  above  equation  into  equation  (7)  gives 
the  following  result. 


f 

osc 


9.269  X 


(11) 


In  order  to  demonstrate  the  sensitivity  of  this  equation  to  the 
various  input  parameters,  we  postulate  a  nominal  pulsed  HEL  system. 
To  construct  such  a  system,  we  assume  that  it  will  be  E-beam  or 
E -beam-controlled  discharge  pumped,  and  operate  at  a  total  pressure 
of  one  atmosphere.  We  further  assume  that  the  active  molecules  are 
one  percent  of  total  and  that  one  percent  of  the  active  molecules  are 
inverted.  Let  us  further  assume  that  we  have  a  duty  cycle  of  0.001 
percent  (i.  e. ,  10  pps,  1  Msec  pulse  width).  With  these  assumptions, 
the  stimulated  emission  rate  can  be  estimated  by  stipulating  that  the 
average  power  must  be  competitive  with  the  DF  HEL.  Thus  the 
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3 

average  output  power  must  be  equivalent  to  the  0.01  kw/cm  of  DF. 
With  those  assumptions  we  obtain: 


S  ^ 


1.879  X  10 
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V 


0 


(12) 


Using  the  nominal  values  other  parameters  of  Table  III  we  obtain  an 

equation  for  f  in  our  nominal  pulsed  HEL  system. 

OS  c 


=  6-75  X  10^^ 
OSC  rr.7/2  2 

^  0 


(13) 


The  resulting  f^^^  values  as  a  function  of  wavenumber  and  tempera¬ 
ture  are  shown  in  Figure  3.  That  figure  rather  clearly  shows  that 
the  shorter  wavelength  systems  are  much  more  susceptible  to  possible 
hole  burning  and  subsequent  gain  saturation  than  are  the  infrared  laser 
systems. 


3.  2  Application  to  Specific  Systems 

In  the  current  period,  we  have  investigated  a  number  of  other 
laser  systems  which  are  of  current  interest.  Equation  11  can  be  used 
to  predict  the  f^^^  directly  for  a  number  of  laser  systems.  The 
principal  difficulty  with  this  equation  is  obtaining  a  reasonable  value 
of  the  stimulated  emission  rate  per  excited  molecule.  Small  scale 
experiments  can  be  used  to  estimate  this  rate  when  such  results  are 
available.  Such  estimates  were  used  in  the  preceding  report  to  derive 
fosc  values  for  a  number  of  promising  HEL  systems.  In  the  present 
period,  we  have  considered  a  broad  range  of  laser  systems.  The 
parameters  used  and  the  resulting  f^^^  values  are  shown  in  Table  III. 
In  determining  the  pulsed  stimulated  emission  rate,  a  one  pulse  per 
second  repetition  rate  at  the  indicated  pulse  width  (w)  has  been 
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Table  III 

Evaluation  of  the  Homogeneous  Broadening  Parameter  f^^^  for  Various  HEL  Systems 


although  such  a  system  may  not  currently  exist. 


assumed.  The  individual  systems  and  the  detailed  results  are  dis¬ 
cussed  in  the  following  paragraphs. 

K„  and  Na^  are  both  members  of  a  class  of  proposed  alkali- 
metal  dimer  laser  systems.  Such  systems  have  been  proposed  by 
York  and  Gallagher  [13]  as  efficient  E-beam  or  discharged  pumped 
systems.  The  K2  system  parameters  given  in  Table  m,  are  similar 
to  those  of  a  system  where  gain,  although  not  lasing,  has  been  mea¬ 
sured  by  Palmer  [14]  in  a  self- sustained  discharge  laser  configuration. 
As  may  be  seen,  the  f  parameter  is  fairly  large  for  this  system 
and  thus  gain  saturation  does  not  appear  to  be  a  problem.  This  result 
should  be  typical  of  the  high  pressure  (1  to  10  atmospheres)  alkali- 
metal  dimer  systems. 

The  parameters  chosen  for  the  Na2  system  are  representative 
of  a  low  pressure  alkali-dimer  system.  This  type  of  laser  is  probably 
less  adaptable  to  high  power  operation  than  the  high  pressure  system, 
however  Na2  has  been  lased  in  a  low  pressure  system  with  optical 
pumping  [15].  It  may  be  seen  that  when  we  assume  HEL  pulsed  sys¬ 
tem  operation  the  f  value  approaches  one.  This  indicates  that 
some  hole  burning  may  take  place  in  such  a  system.  Whether  this  is 
an  important  effect  will  depend  on  the  detailed  operation  of  a  specific 
laser. 

The  third  laser  system  considered  was  a  KrF  laser.  This  laser 
is  one  of  a  class  of  recently  demonstrated  rare-gas  mono-halide  sys¬ 
tems.  In  these  systems  an  E-beam  or  E-beam  sustained  discharge 
is  used  to  excite  the  rare  gas  atom.  This  excited  rare  gas  atom  be¬ 
haves  very  much  like  an  alkali-metal  since  it  now  has  one  electron 
in  its  outer  shell.  Thus  the  alkali-like  rare-gas  reacts  readily  with 
a  halide  to  form  a  strongly  bound  ionic  molecule  in  an  excited  state. 
The  lasing  transition  is  from  the  strongly  bound  excited  state  to  a 

16 


weakly-bound  or  dissociative-ground  state.  Thus  an  inversion  is 
easily  produced.  It  may  be  seen  from  the  table,  that  with  our  as¬ 
sumed  parameters  the  pulsed  high  energy  KrF  laser  system  may 
be  susceptible  to  gain  saturation.  In  fact,  since  it  appears  that  the 
laser  may  have  a  much  shorter  pulse  width  than  has  been  assumed 
in  our  calculation  [16];  inhomogeneous  broadening  may  be  a  serious 
problem  for  a  high  power  laser.  However,  the  low  power  laser  des¬ 
cribed  in  reference  16  does  not  show  this.  If  inhomogeneous  broad¬ 
ening  were  to  occur  for  a  high  power  mono-halide  laser,  it  might  be 
a  serious  problem  since  the  short  (20  nsec)  spontaneous  emission 
lifetime  of  the  upper  state  provides  a  fast  deactivation  channel. 

The  final  system  considered  in  the  present  period  is  the  proposed 
HCl  non-reacting  mixing  laser  [1].  In  that  laser,  the  HCl  would  be 
pumped  by  vibrational  energy  transfer  from  high  temperature  Hg. 

The  resulting  f^g^  parameter  for  this  laser  clearly  shows  that  the 
laser  line  may  be  considered  totally  homogeneous  here.  Thus ,  this 
laser  is  similar  to  the  other  infrared  lasers  considered  in  the  earlier 
report,  where  the  cross -relaxation  rate  is  more  than  adequate  to  com¬ 
pensate  for  the  fact  that  the  homogeneous  Lorentz  widths  are  less  than 
the  inhomogeneous  Doppler  widths. 

In  the  previous  report  a  number  of  other  potential  HEL  systems 
were  also  considered.  The  parameters  used  and  the  resulting  f^^^ 
values  are  given  in  Table  IV.  The  complete  table  has  been  reproduced 
here,  since  a  number  of  minor  errors  occurred  in  the  previous  publica¬ 
tion.  The  final  results,  however  are  essentially  the  same  as  previously 
reported  and  the  conclusions  of  that  report  are  unchanged. 

3. 3  Conclusion 

In  this  period,  we  have  used  the  results  of  a  modified  rate  equa¬ 
tion  theory,  to  assess  the  importance  of  inhomogeneous  broadening  and 
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Table  IV 

Previously  Evaluated  Homogeneous  Broadening  Parameters,  for  Six  Potential  HEL  Systems*  [10] 


o 

o 

o 

CO 

I> 

« 

O 

O 

O 

o 

o 

ID 

O 

O 

X 

o 

ID 

X 

ID 

t- 

LD 

O 

O 

o 

X 

o 

ID 

X 

ID 

tr- 

• 

• 

ID 

O 

CO 

ID 

• 

tH 

C35 

• 

o 

05 

d 

CQ 

CO 

CO 

t-H 

CO 

iH 

t— 1 

o 

ci 

• 

o 

t-H 

CM 

• 

U 

< 

d 


<D 

X! 


Q 


o 

u 


CQ 

<D 

s 

d 

iu 

d 

pL4 


o  o 


CO 
00  I 

o  o 


uo 

o  o 


00 

o  o 


O 

• 

o 

• 

CO 

o 

ID 

• 

O 

O 

ID 

CO 

o 

o 

o 

CO 

CO 

CO 

CM 

'rH 

o 

00 

o 


X 

lO 

tH 


X 

_  o 
•  g  00 

O  O  o 

ir-i  ^ 

CM 


X 

t- 


o  o 


00  o 
o  o 


o 

o 

1> 

05 

• 

O 

o 

o 

CO 

o 

• 

• 

CM 

00 

CM 

CO 

o 

o 

CO 

CO 

iH 

CM 

CM 

LD 

• 

XJ 

o 

CO 

m 


in  § 

o  ^ 


o 

o 


X  X 

O  lO 
CO  O 
m  th 

•  • 


o 

o 

vH 

tH 

ID 

• 

05 

O 

• 

• 

CO 

CO 

ID 

05 

o 

CO 

CO 

tH 

rH 

CO 

CM 

ID 

I  m  CM 

o  o  o 

iH  tH  tH 

XXX 

O  O 

CD  CO 

CM  tH  CO 

•  •  • 

ID  1-H  tH 


CO 

o  o 


CO 


CO 

O 

ID 

O 

X 

o 

X 

tH 

CO 

O 

O 

O 

tH 

LD 

• 

CO 

O 

m 

tr- 

o 

CM 

^  CO 

CM 

CO 

O 

• 

05 

eo 

ID  tr- 

o  o 


o 

rH  CO 

o  o 


X  X 
o  o 

CM  O 
CO 

•  • 


05  CO 

o  o 


o  o 

•  • 

CO  CO 


o<e} 


o<Cj 

h 


cu 

> 

o 

< 


00 

CM 


c 

o 

S3 


5 


o 


X 

X 

X 

X 

o 

CM 

o 

CO 

o 

ID 

CJ5 

ID 

CO 

05 

o 

LD 

• 

O 

CO 

ID 

• 

O 

o 

ID 

CM 

• 

t- 

• 

CO 

CM 

tr- 

tH 

s  s 


S  S' 

H 


§ 

o 

? 

Ph 


C! 

o 

o 


o 

(D 

CO 


T3 

CD 

CO 

4^  S3 

C5  a 


CD 

CQ 

S3 

Oi 


o 

CD  I 


0 

CQ 

i-H 


CO  CM 


•s 

O 

(Q 

CO 

o 

0 

o 

CO 

O 

1 

03 

o 

Xi 

o 

■a 

I 

4-> 

ccJ 


CO 

s 

0) 

4-> 

CO 

>. 

CO 

w 

K 

ol 

•I-l 

4-4 

P3 

CD 

4-> 

o 

o* 


CD 

"c^ 

4-> 

d 

CD 

CQ 

0) 

S4 

a 

CD 

CD 

rO 


CD 

d 

CD 


CO 

CD 

a 

•r-l 

4-4 

CO 

0 

0 

d 

CQ 

0 

4-» 

0 

a 

ct{ 

u 

d 

a 

0 

H 


18 


systems  may  not  currently  exist. 


cross  relaxation  on  the  operation  of  high  energy  laser  systems.  As 
a  result  of  this  work  it  is  clear  that  typical  HEL’s  which  operate  in 
the  infrared  have  a  sufficiently  fast  cross -relaxation  rate  to  insure 
that  the  commonly  used  homogeneous  line  assumptions  are  valid.  On 
the  other  hand,  it  appears  that  the  inhomogeneous  broadening  of  short 
wavelength  pulsed  laser  systems  may  cause  performance  limitation 
under  certain  cavity  conditions.  Thus  if  high  power  visible  lasers 
are  constructed,  the  cavity  conditions  should  be  tailored  to  insure 
sufficient  cross -relaxation  rate  so  that  the  maximum  energy  may  be 
efficiently  extracted  from  the  cavity. 
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4.  0  ALKALI-NOBLE  GAS  LASER  INVESTIGATION 

In  the  present  report,  the  NaAr  molecule  has  been  chosen  as 
prototypical  of  the  alkali-noble  laser  molecules.  Here  we  present 
the  results  of  an  extensive  ab  initio  and  semi -empirical  study  of  the 
potential  functions,  transition  moment  and  spectral  gain  of  this  molec¬ 
ular  system.  The  present  results  indicate  that  the  bound  A  ^  excited 
state  results  mainly  from  an  ion-induced  dipole  interaction  but  with 

possible  contribution  from  Van  der  Waals  forces.  The  transition 

2  2 

moment  has  also  been  calculated  for  the  A  tt  to  X  S  laser  transi¬ 
tion.  This  calculation  indicates  that  the  transition  moment  variation 
with  internuclear  distance  has  a  significant  effect  on  the  gain  of  the 
medium.  Finally  gain  calculations  are  given  for  this  system  and  ex¬ 
trapolation  to  other  Na-noble  gas  systems  is  discussed. 

4. 1  Potential  Curve  Calculations 

The  important  potential  curves  for  the  alkali-noble  laser  systems 

9  2 

are  the  repulsive  X  L  ground  state  and  the  weakly  bound  A  n  first 
excited  state.  Empirical  scattering  measurements  have  given  approxi¬ 
mate  constants  [17]  and  shapes  for  the  ground  states  of  these  systems; 
however  only  recently  have  excited  state  potential  curves  become  avail¬ 
able  [18].  A  series  of  semi-empirical  calculations  have  been  published 
by  Baylis  [19]  and  Pascale  and  Vanderplanque  [20]  for  the  alkali-noble 
system.  However  these  have  consistently  predicted  excited  state  well 
depths  less  than  half  the  values  determined  experimentally. 

2 

4.  2  Upper  Laser  Level  A  n  State 

The  bound  nature  of  alkali -noble  excimers  have  long  been  assumed 
to  be  due  principally  to  the  attractive  Van  der  Waals  forces  which  are 
known  to  be  important  in  long  range  ground  state  interactions  between 
alkali  and  noble  gas  atoms.  For  example,  the  Baylis  theory  which  does 
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explain  greater  than  fifty  percent  of  the  well  depth  for  most  alkali-noble 

systems  is  based  on  the  extrapolation  of  the  Van  der  Waals  attractive 

forces  to  the  excited  molecular  states.  Thus,  in  order  to  calculate 

more  accurately  the  configuration  interaction  energies,  including 

2 

Van  der  Waals’  contribution  of  the  A  v  state  of  NaAr,  an  extensive 
ab  initio  calculation  was  carried  out. 

2 

4.  2. 1  A  IT  Configuration  Interaction  Results 

In  the  NaAr  Configuration  Interaction  (Cl)  calculation,  we  have 
followed  a  MC-SCF  variational  method  previously  used  by  Hirchfelder 
[21]  and  Das  [22].  In  this  approach,  a  nominal  size  basis  set  is  used 
to  determine  the  configuration  interaction  energies  between  the  prin¬ 
cipal  Hartree-Fock  (HF)  configuration  and  the  dispersion  and  other 
configurations  which  are  expected  to  give  significant  attractive  force 
contributions.  It  should  be  emphasized  however,  that  only  the  Cl  po¬ 
tential  contributions  have  been  calculated  here,  since  the  basis  set 
used  is  not  large  enough  to  predict  the  total  interaction  energy  of  the 
state  accurately.  The  basis  set  and  configurations  used  are  listed  in 
Tables  V  and  VI.  From  Table  V,  it  can  be  seen  that  a  number  of 
configurations  which  were  omitted  from  our  previous  report  [10]  be¬ 
cause  of  certain  limitations  of  the  computer  codes  have  been  included 
here  through  use  of  an  alternate  program,  ♦  The  results  of  this  cal¬ 
culation  are  given  in  Table  VII.  The  3pj^^  -  3s^^  results  are  omitted 
here,  because  in  the  MC-SCF  itteration  these  configurations  became 
identical  with  configurations  contained  in  the  3pj^^  -  3p^  group. 

It  is  useful  to  analyse  these  results  in  terms  of  the  expected 
Cg/R  dependence  of  the  dispersion  energy  and  the  expected  C^/k 

*  We  have  recently  been  informed  that  these  limitations  have 

since  been  removed. 
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Table  V. 

2 

Dominant  Configurations  Used  in  the  Calculation  of  the  A  n  State  of  NaAr 


Description 

Configurations 

No.  of 

Vector  Couplings 

Hartree-Fock 

1 

(5PNa>  ^PAr> 

^®’^Na  2  2  3 

3dn^^(3sa'^3p<T  3piT  3d7r)^ 

* 

4,  2 

3d6Na(3s(T^3pa^3p;r^3dCT)Ar 

2"^ 

3daNa(3s<^^3p^^3Pff  ^Sda)  ^ 

2 

3dCTNa(3sa^3pCT^3p77^3d6)Ar 

2'*' 

3d6Na(3st^^3pa^3p7r^3d6)Ar 

2,  2" 

* 

3d6  Na(2 

2“^ 

2  4 

3daNa(3sa  3pCT3d7T3pi7 

2 

3d7rNa(3sa^3pa3da3p7r^)Ar  . 

2 

<3PNa> 

2  4 

3d7Tj^^(3sa3po-  4pa3pff 

2 

3d6^^(3sa3pa^4p7r3pir^)  ^ 

2"^ 

3dCTj^J3sa3pa^4pjr3pir^)^ 

2 

Overlap-Charge- 

7 

Transfer 

7 

*  These  omitted  configurations  will  contribute  only  to  dispersive  terms 
of  higher  order  than  the  Cg  term. 

+  The  contributions  from  these  configurations  were  calculated  using 

MC-SCF  codes  provided  by  Dr.  J.  Hinze  of  the  University  of  Chicago. 
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Table  VI.  2 

Basis  Set  Used  for  the  Calculation  of  the  A  n  State  of  NaAr 


nl 

a -Set 
Exponent 

nl 

TT-Set 

Exponent 

nl 

6 -Set 
Exponent 

Na  10 

10.6259^ 

Na  21 

3.4009 

Na  32 

0.75 

20 

3.  2857^ 

31 

0.8592 

42 

0.75 

30 

1.24^ 

31 

0. 5292 

43 

0.8592 

30 

0.  75^ 

32 

0.75 

43 

0.  5292 

21 

3.  4009^ 

42 

0.  75 

31 

0.8592^ 

43 

0. 8592 

31 

0. 5292^ 

43 

0.  5292 

32 

0.75^ 

42 

0.75^^ 

1 

43 

0.  8592*^ 

43 

0.  5292^ 

Ar  10 

17.  5057^ 

At  21 

7.0041 

Ar  32 

2.  893 

20 

6. 1152^ 

31 

2.893 

32 

1.605 

30 

2.  5856^ 

31 

1.605 

42 

2.893 

21 

7.0041^ 

41 

2.  5856 

42 

1.605 

31 

2.  893^^ 

32 

2.  893 

43 

2.  2547 

31 

1.  605^ 

32 

1.605 

41 

2.  5856^^ 

42 

2.  893 

32 

2.893*^ 

42 

1.605 

32 

1.  605*^ 

43 

2.  2547 

42 

2.893^ 

42 

1.605^ 

43 

2.  2547® 

a.  Taken  from  either  the  minimal  or  nominal  basis  set  for  Na  or  Ar  as  listed 
in  the  BISON  program  of  Argonne  National  Laboratory. 

b.  Taken  from  the  HF  calculation  for  Na-Ar  by  Janis  and  Wahl  [25]. 

c.  Chosen  for  the  excitations  as  prescribed  in  Ref.  22. 
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Table  VII. 

Results  of  Configuration  Interaction  Calculation  for  the 

State  of  NaAr 


Internuclear 

Separation 

(Bohr) 

^PNa-^PAr 
(Hartree  x  10^) 

Charge  - 
Transfer 

(Hartree  x  10^) 

^CI 

(Hartree) 

5.5 

-1.  966 

-0.  265 

-2.321 

6.0 

-1.604 

-0. 198 

-1.802 

7.0 

-1.005 

-0. 114 

-1. 119 

10.0 

-0.224 

-0.  012 

-0.236 
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dependence  for  the  ion-induced  dipole  potential.  The  Cg  coefficient 
determined  from  the  Cl  results  is  shown  as  a  function  of  R  in  Fig¬ 
ure  4.  The  ten  bohr  Cg(R)  value  of  236  in  atomic  units  is  slightly 
larger  than  the  calculated  asymptotic  value  of  190  for  the  ground 
state  of  NaAr  [23],  The  expected  decrease  of  the  Cg(R)  coefficients 
as  internuclear  distance  decreases  is  similar  to  the  results  of  pre¬ 
vious  calculations  on  similar  systems  [24]. 

The  Cl  energy  contribution  has  also  been  used  to  determine  the 
C^(R)  coefficient  expected  with  an  ion-induced  dipole  interaction.  The 
results  of  such  a  calculation  are  given  in  Figure  4.  From  that  figure 
it  is  clear  that  the  C^(R)  analysis  gives  a  much  flatter  R  dependence, 
leading  to  the  conclusion  that  the  principal  attractive  forces  here  are 
not  dispersion  forces  but  rather  the  ion-induced  dipole  attraction. 

2 

4.  2.  2  An  Hartree-Fock  Potential 

As  discussed  in  our  previous  report,  it  was  expected  that  the  ad¬ 
dition  of  the  calculated  Cl  potential  to  the  large  basis  HF  potential 

computed  by  Janis  and  Wahl  [25]  would  provide  an  accurate  represen- 

2 

tation  of  the  NaAr  A  n  potential.  However  after  a  thorough  examina¬ 
tion  of  results  obtained  for  the  small  basis  HF  configuration  included 
in  the  current  MC-SCF  calculation,  it  became  apparent  that  an  ion- 
induced  dipole  interaction  is  also  important  for  the  HF  potential.  In 

fact,  in  contrast  to  the  earlier  HF  results,  our  HF  results  show  sig- 

2  2 

nificant  binding  due  to  a  mixture  between  the  A  n  and  higher  n  states. 
The  dominant  effect  seen  is  the  polarization  of  the  Na  3piT  orbital  away 
from  the  incoming  Ar  atom  through  orbital  mixing  with  the  Na  4d'Tr 
orbital.  Mixing  is  also  observed  with  the  3pir  and  4p‘n-  Ar  orbitals. 
Thus  the  valence  orbital  is  taking  on  the  character  of  a  united  atom 
orbital  and  the  Na  atom  begins  to  take  on  the  character  of  the  Na^ 
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ion  from  the  perspective  of  the  Ar  atom.  The  previous  HF  calcula¬ 
tion  did  not  show  this  binding  effect  because  neither  the  important 
4pTr  nor  the  4d'Tr  orbitals  were  included  in  the  basis  set  used  there. 

On  the  other  hand,  the  basis  set  used  in  the  present  MC-SCF  calcula¬ 
tion  is  simply  too  small  to  give  an  accurate  HF  potential.  Neither  of 

these  HF  potentials  now  appears  to  be  an  accurate  representation  of 
2 

the  true  A  it  state  HF  potential  in  light  of  our  present  understanding 
of  the  physics,  but  the  large  basis  HF  potential  can  conveniently  be 
used  with  the  experimental  results  to  get  a  better  understanding  of 
the  importance  of  the  attractive  potential  effects. 

2 

4.  2.  3  NaAr  A  it  -  Total  Attractive  Potential 

The  attractive  potential  for  intermediate  internuclear  separations 
is  conveniently  obtained  by  substracting  the  HF  energy  from  the  ex¬ 
perimental  potential  values  of  York  [18].  When  this  is  done  and  the  re- 
suits  are  analysed  as  or  Cg/R°  interactions,  the  C(R)  curves  of 

Figure  5  are  obtained.  Again,  from  that  figure  it  is  apparent  that  the 
interpretation  gives  the  more  constant  value  of  the  force  constant. 

It  seems  clear  from  our  Cl  results  and  the  semi-empirical  analysis 
given  above,  that  in  the  well  region  near  the  potential  minimum  the  prin¬ 
cipal  bonding  force  results  from  an  ion-induced  dipole  effects.  The  turn 
over  of  the  C^(R)  in  Figure  4,  beyond  8  bohrs  indicates  that  at  that  point 
the  Cg  dispersion  forces  probably  are  becoming  important  while  the 

ion-induced  dipole  forces  are  becoming  less  significant.  The  calculated 

2 

Cl  contribution  and  the  total  potential  of  the  A  ft  state  is  given  in  Fig¬ 
ure  6  and  Table  VIII. 


4.  3  Lower  Laser  Level  -  NaAr  X 
2 

The  X  potential  is  important  because  the  shape  and  relative 

position  of  this  potential  with  respect  to  the  excited  state  determines 
the  peak  laser  gain  wavelength.  Two  types  of  measurements  are 
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Ce(R)  (Hartree 


C^(R)  (Hartree-Bohr 


Table  Vin.  Various  Contributions  to  the  A  n  Potential  of  NaAr 


R(AU) 

Vhf[25] 

'^CI 

* 

'^Att 

V 

^  total 

V  [18] 
exp’^  ■' 

Bohr 

Hartree  x  10 

3 

Hartree  x  10 

Hartree  x  10 

Hartree  x  10^ 

-1 

cm 

-1 

cm 

5.0 

4.  35 

-2.  816 

0.  526 

115 

5.5 

1.69 

-2.  231 

-3.461 

-1.771 

-389 

-389 

6.0 

0.  79 

-1.802 

-3.  229 

-2.439 

-535 

-517 

6.5 

0.  42 

-1.405 

-2.  626 

-2.  206 

-484 

-498 

7.0 

0.  26 

-1. 119 

-2. 046 

-1. 786 

-392 

-395 

7.5 

0. 17 

-0.876 

-1.  579 

-1.409 

-309 

-302 

8.0 

0. 13 

-0.687 

-1.230 

-1. 100 

-241 

-238 

8.5 

0.  07 

-0.  530 

-0.  886 

-194 

9.0 

0.  05 

-0.406 

-0.  695 

-153 

9.5 

0.  02 

-0.  306 

-0.559 

-123 

10.0 

-0.  00 

-0.  236 

-0.  459 

-100 

*  Contains  V^j  contributions. 
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Table  I X.  Various  Contributions  to  the  X  S  Potential  of  NaAr 


R 

Vhf[25] 

^disp 

V 

total 

Bohr 

- 7^ 

Hartree  x  10 

.  _3 

Hartree  x  10 

Hartree  x  10 

-1 

cm 

Hsm 

13.19 

-2.83 

10.36 

2274 

2375 

6.0 

5.86 

-1.72 

4. 14 

909 

795 

2.70 

-1.07 

1.63 

358 

315 

1.09 

-0.  64 

0.  45 

99 

142 

0.  40 

-0.  35 

0.05 

11 

0.  07 

-0. 19 

-0. 12 

-26 

0 
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available  for  this  potential,  but  both  have  certain  limitations.  Scat¬ 
tering  measurements  results  of  Malerick  and  Cross  [17]  are  inter- 

fi 

preted  as  giving  a  Cg/R  repulsive  potential  with  a  value  of  the  re¬ 
pulsive  constant  of  3.71  in  atomic  units.  This  however,  must  be 
treated  as  only  an  approximation  to  the  true  curve.  The  other  experi¬ 
mental  potential  curve  was  determined  by  York  [18]  from  pressure 
dependent  emission  measurements.  The  shape  and  position  of  this 
potential  are  not  unambiguously  determined  by  the  fluorescence  mea¬ 
surement  since  the  analysis  of  those  results  requires  that  a  starting 
radius  to  fix  the  potential  position  be  taken  from  another  measure¬ 
ment. 


4.3. 1 


Hartree-Fock  Potential 


2 

The  Hartree-Fock  potential  for  the  X  ^^^^2  state  of  NaAr  has 

been  calculated  by  Janis  and  Wahl  [25]  with  the  same  basis  set  used 
2 

for  the  A  state.  For  this  state,  the  ion-induced  dipole  effect  should 
be  negligible  because  of  the  symmetry  of  the  3pcr  state.  Thus  we 
would  expect  the  HF  calculation  to  give  an  accurate  representation  of 
the  single  configuration  energy.  Of  course,  since  configuration  inter¬ 
action  is  not  included  here,  dispersion  effects  must  be  added  to  the 
HF  results  to  give  an  accurate  representation  of  the  true  molecular 
potential. 


4.  3.  2  Waals  Potential 

The  Van  der  Waals  calculation  for  this  state  has  not  been  com¬ 
pleted;  however,  art  approximate  Van  der  Waals  potential  for  this  state 
can  be  determined  through  the  use  of  HF  results  and  the  scattering 
measurement.  The  resulting  Cg(R)  coefficient  values  from  such  a 
calculation  are  plotted  in  Figure  7.  We  have  simply  assumed  the 
smooth  curve  indicated  on  that  figure  for  the  dispersion  coefficient 
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and  combined  the  resulting  dispersion  potential  with  the  HF  potential 

2 

of  Janis  to  obtain  the  X  L  potential  given  in  Table  DC  and  Figure  8.  It 
can  be  seen  from  that  figure  that  this  potential  is  in  reasonable  agree¬ 
ment  with  the  potential  obtained  by  York  [19], 

4.4  Transition  Moment  Calculations 

The  dipole  transition  moment  is  also  required  for  laser  gain  cal¬ 
culation.  Earlier  computations  have  assumed  that  the  transition  mo¬ 
ment  is  independent  of  the  internuclear  separation  [10,  26].  In  the 
present  study  we  have  calculated  the  transition  moment  in  the  HF 
approximation  with  the  wave  functions  determined  by  Janis  and  Wahl. 
From  those  wave  functions,  both  the  dipole  length  and  dipole  velocity 
transition  moments  forms  were  calculations  with  the  use  of  an  Argonne 
computer  code  [27].  The  resulting  oscillator  strengths  and  transition 
moment  values  for  NaAr  are  given  in  Table  X.  These  results  may  be 
compared  with  the  dipole  length  oscillator  strength  values  of  Table  XI 
computed  by  Krauss,  Maldonado  and  Wahl  [27]  for  the  similar  NaHe 
molecule.  The  divergence  between  the  f  ^  and  f^  forms  of  the  oscillator 
strength  is  disturbing,  however,  in  light  of  the  similarity  between  the  f ^ 
values  for  NaAr  and  NaHe,  these  f ^  values  might  be  preferred.  The  ion- 
induced  dipole  and  the  Van  der  Waals  wave  functions  would  also  be  ex¬ 
pected  to  provide  significant  correction  to  the  HF  transition  moment, 
so  the  values  given  here  must  be  considered  approximate. 

4.  5  NaAr  Gain  Calculation 

Under  the  quasi-static  theory  assumption,  the  reduced  laser  gain 
is  given  by  the  following  expression  [28]: 

© 


& 


g(v)  = 


u  R 


& 


fl 


1^1 
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Table  X, 


Calculated  Oscillator  Strengths  and  Dipole  Transition  Moments 

O  _L  2 

for  the  X  E  -  A  it  Transition  of  NaAr 


Inter  nuclear 
Separation 
R(bohr) 

Oscillator  Strength 

Square  of  Transition 
Moment  1^2 
(electron  -  bohr)2 

(dipole  length) 

(dipole  velocity) 

f 

V 

(from  fj) 

(from  f^) 

5.0 

0.  578 

0.  666 

6.82 

7.  86 

6.  0 

0.607 

0.  665 

6.76 

7.  40 

8.0 

0.639 

0.  656 

6.71 

6.  87 

OO 

Exp.  f 

=  0.654 

=  6.35 

Table  XI 

Oscillator  Strength  and  Dipole  Transition  Moment 
of  NaHe  [27]  X^5  -  A^ir 


Internuclear 

Separation 

R(bohr) 

Oscillator  Strength 
(dipole  length) 

Square  of  Transition 
Moment 

(electron  -  bohr)'^ 

5.0 

0.610 

6.92 

6.0 

0.631 

6.87 

8.0 

0.  646 

6.73 

OO 

0.  649 

6.72 

*  Calculated  from  the  fj^  given  in  reference  [27]. 
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The  frequency  of  the  transition  is  given  by 
i^(R)  =  V'(R)  -  V'’(R). 

For  the  lower  and  upper  state  potentials  V”  and  V',  lease  square 
fits  to  the  results  given  in  Tables  VIII  and  DC  were  used. 

The  Einstein  coefficient  is  represented  by  the  expression 


A^l(v)  =  A(Xq) 


(R)  2 


where  u  ^(R)  is  fit  to  the  results  of  Table  X  and  is  the  expert- 
mental  oscillator  strength  for  th6  Na  Dg  transition.  A  value  of  A(Xq) 
of  6.667  X  10*^  sec~^  was  used  with  Xq  equal  to  5890  A,  A  translational 
temperature  of  650°  K  and  an  electronic  temperature  of  7000°  K  have 
been  assumed. 


The  resulting  reduced  gain,  stimulated  emission  and  absorption 
coefficients  are  given  in  Figures  9  and  10.  In  Figure  9,  the  constant 
value  of  the  transition  moment  has  been  assumed.  This  leads  to  a  peak 
gain  of  1.24  x  lO"^^  cm^  at  6640  k.  This  wavelength  corresponds  to 
an  inter-nuclear  separation  of  5.55  bohrs  which  is  on  the  inner  wall  of 
the  excited  state  potential.  Figure  10  shows  the  range  of  results  ob¬ 
tained  with  the  two  forms  of  the  calculated,  radius-dependent  transition 
moment  representations.  The  f^  form  gives  the  largest  peak 
gain  value  of  1.69  x  10"^^  cm^  almost  40%  greater  than  the  constant 
moment  result. 


4.  6  Conclusion 

A  number  of  alkali-noble  gas  laser  systems  have  been  investigated 

2 

theoretically.  This  investigation  has  indicated  that  the  A  tr  excited 
bound  state  of  the  NaAr  molecule  is  basically  a  Rydberg  state  where 
ion-induced  dipole  and  Van  der  Waals  forces  both  significantly  contribute 
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Reduced  Gain/Stimulated  Emission/Absorption  (cm  ) 


Wavelength  (^) 


Figure  9.  Reduced  Gain,  Stimulated  Emission  and  Absorption  for 
NaAr  with  a  Constant  Transition  Moment 
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Wavelength  (A) 


Figure  10.  Reduced  Gain,  Stimulated  Emission  and  Absorption  for 
NaAr  Calculated  with  a  R-Dependent  Transition  Moment 


; 


to  the  binding.  In  the  course  of  this  investigation  a  number  of  theo¬ 
retical  gain  calculations  which  can  be  used  to  predict  the  gain  for  actual 
laser  media  conditions  have  been  performed.  It  has  been  shown  that 
predicted  variations  of  the  transition  moment  can  cause  a  significant 
change  in  the  predicted  gain  in  the  NaAr  system. 

The  current  study  indicates  that  ab  initio  molecular  calculations 
are  useful  in  obtaining  information  on  the  chemistry  of  alkali -noble 
molecule  formation  important  for  laser  design.  However,  such  cal¬ 
culations  are  expensive,  and  care  must  be  exercised  to  limit  the  quan¬ 
tity  of  calculations,  to  include  the  minimum  required  to  obtain  the 
desired  theoretical  information.  We  find  that  the  physical  insight 
gained  from  ab  initio  calculations  on  intermediate  size  systems  like 
NaAr  is  considerable,  and  that  it  can  be  applied  with  more  confidence 
to  interpret  larger,  similar  molecules  than  can  extrapolations  from  the 
more  typical,  small  molecule  ab  initio  results. 
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AMCPM-HEL-T  (Dr.  Evers) 

Commander,  US  Army  Missile  Command 
Attn:  AMSMI-RNS 

Commanding  Officer 

US  Army  Mobility  Equipment  R&D  Center 
Attn:  AMEFB-MW 


Commander,  Rock  Island  Arsenal 
Attn:  SARRI-LR 

Mr.  J.  W.  McGarvey 

Commander,  US  Army  Armament  Command 
Attn:  AMSAR-RDT 

Director 

Ballistic  Missile  Defense  Advance 
Technology  Center 

Attn:  ATC-O,  Mr.  W.  O.  Davies 

Commander,  USA  Test  &  Evaluation 
Command 

Attn:  AMSTE-ME  (Dr.  N.  Pentz) 

Commander,  US  Army  Materiel  Command 
Attn:  AMCRD-T  (Mr.  Paul  Chernoff) 

(Dr.  David  Stefanye) 
(Dr.  B.  Zarwyn) 

Director,  US  Army  Ballistic  Research  Lab. 
Attn:  Dr.  Robert  Eichelberger 
Mr.  Frank  Allen 
Dr.  E.  C.  Alcarez 

Commandant,  US  Army  Air  Defense  School 
Attn:  Air  Defense  Agency 

Commandant,  US  Army  Air  Defense  School 
Attn:  ATSA-CD-MS 

Commander,  US  Army  Training  and 
Doctrine  Command 
Attn:  ATCD-CF 

Commander,  USA  Frankford  Arsenal 
Attn:  Mr.  M.  Elnick  (SARFA-FCD) 

Commander,  US  Army  Electronics  Comman 
Attn:  AMSEL-CT-L  (Dr.  R.  G.  Buser) 

Commander 

US  Army  Combined  Arms  Combat 
Developments  Activity 

Deputy  Commandant  for  Combat  & 

Training  Developments 
US  Army  Ordnance  Center  and  School 
Attn:  ATSL-CTD-MS-R  (LTC  Stewart) 

Department  of  the  Navy 
Office  of  the  Chief  of  Naval  Operations 
Attn:  CDR  L.  E.  Pellock,  USN 
(OP-982F3) 

Mr.  L.  E.  Triggs  (OP-35E) 
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Office  of  Naval  Research,  Boston 
Attn:  Dr,  Fred  Quelle 

Office  of  Naval  Research,  Virginia 
Attn:  Dr.  W.  J.  Condell  (421) 

Department  of  the  Navy 

Deputy  Chief  of  Naval  Materiel  (Dev) 

Attn:  Mr*  R*  Gaylord  (MAT  032B) 

Naval  Missile  Center 

Attn:  Gary  Gibbs  (Code  5352) 

Commander,  Naval  Sea  Systems  Command 

Department  of  the  Navy 

Attn:  Capt.  J.  G.  Wilson,  PMS-405 

Superintendent,  Naval  Postgraduate  School 
Attn:  Library  (Code  2124) 

tJS  Naval  Weapons  Center,  China  LaLe 
Attn:  Mr.  E.  B.  Niccum  (Code  5114) 

Naval  Research  Lab. ,  Washington,  D.  C. 

Attn:  Dr.  J.  M.  MacCallum  (Code  5503) 

EOTPO 

Dr.  P.  Livingston  (Code  5560) 

Mr.  D.  J.  McLaughlin  (Code  5560) 

Dr.  J.  L,  Walsh  (Code  5503) 

Dr.  J.  T.  Schriempf  (Code  6410) 

Dr.  R.  F.  Wenzel  (Code  6410) 

‘  Mr.  R.  W.  Rice  (Code  6360) 

Dr.  L.  R.  Hettche  (Code  6310) 

Dr.  J.  K.  Hancock  (Code  6110) 

Naval  Surface  Weapons  Center,  White  Oak  Lab. 
Attn:  Dr.  E.  L.  Harris  (Code  313) 

Dr.  L.  H.  Schindel  (Code  310) 

Mr.  D.  L.  Merritt  (Code  034) 

Mr.  J.  Wack  (Code  048) 

Hq.  AFSC/XRLW,  Andrews  AFB 
Attn:  Major  James  M.  Walton 

Hq.  USAF(RDPS),  Washington,  D.  C. 

Attn:  LTCA.  J.  Chiota 

Hq.  AFSC(DLCAW),  Andrews  AFB 
Attn:  Maj.  H.  Axelrod 

Air  Force  Weapons  Lab. ,  Kirtland  AFB 
Attn:  Col.  Donald  L.  Lamberson  (AR) 

Col.  John  C.  Scholtz  (PG) 

Col.  Russell  K.  Parson  (LR) 

Lt.  Col.  John  C.  Rich  (AL) 


Hq.  SAMSO,  Los  Angeles 

Attn:  Capt.  Dorian  A*  DeMaio  (XRTD) 
IND 

Capt.  Thomas  J.  Ernst  (DYD) 

AF  Avionics  Lab.  (TEO),  Wright 
Patterson  AFB 

Attn:  Mr.  K.  Hutchinson 

AF  Materials  Lab. ,  Wright  Patterson  AFB 
Attn:  Maj.  Paul  Elder  (LPJ) 

Dr.  William  Frederick  (LPO) 

Department  of  the  Air  Force 
Hq.  Foreign  Technology  Div. 

Attn:  Mr.  K  W.  Buxton  (ETEO) 

AF  Aero  Propulsion  Laboratory 

Attn:  Maj.  George  Uhlig  (AFAPL/NA) 

RADC  (OCSE/Mr.  R,  Urtz),  Griffiss  AFB 

Hq,  Electronics  Systems  Div.  (ESD) 

Attn:  Capt.  Allen  R.  Tobin  (XRE) 

AF  Rocket  Propulsion  Lab. ,  Edwards  AFB 
Attn:  B.  R,  Bornhorst  (LKCG) 

CINCSAC/INEP,  QffuttAFB 

USAF/INAKA,  Washington,  D.  C. 

Attn:  LTC  Frederic  C.  Dunlap 

Defense  Intelligence  Agency,  Washington, 
Washington,  D.  C. 

Attn:  Mr.  Seymour  Berler  (DTIA) 

Central  Intelligence  Agency, 

Washington,  D.  C. 

Attn:  Mr.  Julian  C.  Nall  (OSI/PSTD) 

Dr.  John  E.  Ashman  (OWI/DSD) 

Aerodyne  Res.  Inc. ,  Burlington,  MA 
Attn:  Charles  E.  Kolb 

Analytic  Services,  Inc. ,  Virginia 
Attn:  Dr.  John  Davis 

Aerospace  Corp. ,  Los  Angeles 
Attn:  Dr.  G.  P.  Millburn 

Dr.  Walter  PL  Warren,  Jr. 

Dr.  Elliott  L.  Katz 

Airesearch  Manuf.  Co. ,  Los  Angeles 
Attn:  Mr,  A.  Colin  Stancliffe 
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Johns  Hopkins  University 
Attn:  Dr.  Albert  M.  Stone 
Dr.  R-  E.  Gorozdos 


Atlantic  Res.  Corp. ,  Virginia 
Attn:  Mr.  Robert  Naismith 

AVCO  -  Everett  Res.  Lab. ,  MA 
Attn:  Dr.  George  Sutton 
Dr.  Jack  Dougherty 

Battelle  Columbus  Laboratories 
Attn:  Mr.  Fred  Tietzel  (STOIAC) 

Bell  Aerospace  Co. ,  Buffalo,  NY 
Attn:  Dr.  Wayne  C.  Solomon 

Boeing  Co. 

Attn:  Mr.  M.  I.  Gamble 

ORGN  2-5006,  MS  8C-88 

ELS  Inc. ,  Sunnyvale,  CA 

Attn:  Mr.  Harold  A.  Malliot 

Electro-Optical  Systems,  Pasadena,  CA 
Attn:  Dr.  Andrew  Jensen 

General  Electric  Co. ,  Pittsfield,  MA 
Attn:  Mr.  D.  G.  Harrington 
Room  1044 

General  Research  Corp. ,  Santa  Barbara,  CA 
Attn:  Dr.  R.  Holbrook 

General  Research  Corp. ,  VA 
Attn:  Dr.  Giles,  F.  Crimi 

Hercules,  Inc. ,  Wilmington,  DE 
Attn:  Dr.  R.  S.  Voris 

Hercules,  Inc. ,  Cumberland,  MD 
Attn:  Dr.  Ralph  F.  Preckel 

Hughes  Research  Labs. ,  Malibu,  CA 
Attn:  Dr.  D.  Forster 

Dr.  Arthur  N.  Chester 
Dr.  Viktor  N.  Chester 
Dr.  Gerald  S.  Picus 

Hughes  Aircraft  Co. ,  Culver  City,  CA 
Attn:  Dr.  Eugene  Peressini  (Bldg.  6, 
MS/E-125) 

Dr.  John  Fitts  (MS  5B-138) 

Dr.  J.  A.  Alcalay  (Bldg.  6, 
MS/E-182) 

Hughes  Aircraft  Co. ,  Fullerton,  CA 
Attn:  Dr.  William  Yates 

Institute  for  Defense  Analyses,  VA 
Attn:  Dr.  Alvin  Schnitzler 


Lawrence  Livermore  Lab. ,  Livermore,  CA 
Attn:  Dr.  R.  E.  Kidder 
Dr.  E.  TeUer 
Dr.  Joe  Fleck 
Dr.  John  Emmett 
Mr.  Carl  Haussmann 

Los  Alamos  Scientific  Labs. 

Attn:  Dr.  Keith  Boyer  (MS  530) 

Dr.  O.  P.  Judd 

Lulejian  &  Associates,  Inc. ,  Torrance,  CA 

Lockheed  Palo  Alto  Res.  Lab. 

Attn:  L.  R.  Lunsford 

Orgn.  52-24,  Bldg.  201 

Mathematical  Sciences  Northwest,  Ihc. , 

BeUevue,  WA 

Attn:  Mr.  Peter  H.  Rose  & 

Mr.  Abraham  Hertzberg 

Martin  Marietta  Aerospace,  Denver,  CO 
Attn:  Mr.  Roy  J.  Heyman  (Mail  No.  048 
Dr.  Scott  GiUes  (Mail  No.  8105) 

Massachusetts  Inst,  of  Technology 

Lincoln  Lab. 

Attn:  Dr.  S.  Edelberg 

Dr.  L.  C.  Marquet 
Dr.  J.  Freedman 
Dr.  G.  P.  Dinneen 
Dr.  R.  H.  Rediker 

McDonnell  Douglas  Astronautics  Co. 

Himtington  Beach,  CA 

Attn:  Mr.  P.  L.  Klevatt 

.  Dept.  A3-360-B3GO,  M/S  14-1 

McDonnell  Douglas  Res.  Labs, 

St.  Louis,  MO 

Attn:  Dr.  D.  P.  Ames 

MITRE  Corp. ,  Bedford,  MA 
Attn:  Mr.  A.  C.  Cron 

Dr.  Anthony  N.  Pirri 

Physical  Sciences  Inc. ,  Wakefield,  MA 

Northrop  Corp. ,  Hawthorne,  CA 
Attn:  Dr.  Gerald  Hasserjian 
Dr.  M.  M.  Mann 
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pacific  Sierra  Res.  Corp. 

Attn:  Dr.  TL  Lutomirski 

Philco  Ford,  Newport  Beach,  CA 
Attn:  W.  H.  Rohrer 

Air  Strike  Programs 

RAND  Corp. ,  Santa  Monica,  CA 
Attn:  Dr.  Claude  R  Culp 

Raytheon  Co. ,  Waltham,  MA 

Attn:  Dr.  Frank  A.  Horrigan  (Res.  Div. ) 

Ra}i:heon  Co. 

Bedford  Labs. ,  Missile  Systems  Div. 

Attn:  Dr.  H.  A.  Mehlhorn 
Optical  Systems  Dept. 

M/S  S4-55 

Raytheon  Co. ,  Waltham,  MA 
Attn:  Dr.  Hermann  Statz 

Radio  Corporation  of  America 
Missile  and  Surface  Radar  Div. 

Attn:  Mr.  J.  A.  Colli gan 
Information  Control 

Riverside  Research  Institute,  NY 
Attn:  Dr.  L.  H.  Neill 
Dr.  John  Bose 
HPEGL  Library 

R&D  Associates,  Inc. ,  Santa  Monica,  CA 
Attn:  Dr.  R.  E.  LeLevier 
Dr.  R  Hundley 

Rockwell  International  Corp, ,  Anaheim,  CA 
Attn:  R  E.  Hovda  (DB29) 

Dr.  J.  Winocur  (D/528,  HA14) 

Rockwell  International  Corp. 

Rocketd5me  Div. 

Albuquerque,  NM 

Attn:  Mr.  C.  K.  Kraus,  Manager 

SANDIA  Labs.  Albuquerque,  NM 
Attn:  Dr.  A.  Narath,  ORG  5000 

W.  J.  Schafer  Associates,  Inc. 

Wakefield,  MA 

Attn:  Francis  W.  French 

Stanford  Research  Institute 
Menlo  Park,  CA 

Attn:  Dr.  R  A.  Armistead 
Mr.  J.  E.  Malick 


Science  Applications,  Inc. ,  La  Jolla,  CA 
Attn:  Dr.  John  Asmus 

Science  Applications,  Inc. ,  Arlington,  VA 
Attn:  Mr,  Lawrence  Peckham 

Science  Applications,  Inc. ,  Ann  Arbor,  MI 
Attn:  Dr.  R  E.  Meredith 

Science  Applications,  Inc. ,  Bedford,  MA 
Attn:  Dr.  Robert  Greenberg 

Systems  Consultants,  Inc.,  Washington,  D. C 
Attn:  Dr.  Robert  B,  Keller 

Systems,  Science  &  Software,  La  Jolla,' GA 
Attn:  Mr.  Alan  F.  Klein 

Thiokol  Chemical  Co. ,  Wasatch  Div. 

Attn:  Mr.  James  E.  Hansen 

TRW  Systems  Group,  Redondo  Beach,  CA 
Attn:  Mr.  Norman  F.  Campbell 
Mr.  Eugene  M.  Noneman 

United  Aircraft  Res.  Lab. ,  East  Hartford, 
Attn:  Mr.  G.  H.  McLafferty 
Mr.  Albert  Angelbeck 

United  Aircraft  Corp. ,  Pratt  and  Whitney. 
Acft.  Div. ,  Florida  R&D  Center 
Attn:  Dr.  R  A.  Schmidtke 
Mr.  Ed  Pinsley 

VARIAN  Associates,  EIMAC  Div. 

San  Carlos,  CA 

Attn:  Mr.  Jack  Quinn 

Vought  Systems  Div. ,  LTV  Aerospace  Corp 
Dallax,  TX 

Attn:  Mr.  F.  G.  Simpson 
■  Mail  Station  2-54142 

Westinghouse  Electric  Corp. , 

Defense  and  Space  Center, 

Baltimore,  MD 

Attn:  Mr.  W.  F.  List 

Westinghouse  Res.  Lab. 

Pittsburgh,  PA 

Attn:  Dr.  E.  P.  Riedel 

Mr.  R  L.  Hundstad 
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